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The total lowest-order electromagnetic radiative corrections to the observables in M0ller scatter- 
ing of longitudinally polarized electrons have been calculated. The final expressions obtained by 
the covariant method for the infrared divergency cancellation are free from any unphysical cut-off 
parameters. Since the calculation is carried out within the ultrarelativistic approximation our result 
has a compact form that is convenient for computing. Basing on these expressions the FORTRAN 
code MERA has been developed. Using this code the detailed numerical analysis performed under 
SLAC (E-158) and JLab kinematic conditions has shown that the radiative corrections are significant 
and rather sensitive to the value of the missing mass (inelasticity) cuts. 



PACS numbers: 12.15.Lk, 13.88.+e, 25.30.Bf 



I. INTRODUCTION 



The present intense interest of physicists in polarized 
M0ller scattering is stimulated by several reasons. To- 
day the measurement of the parity-violating asymmetry 
Apv in the recent experiment E158 0,0 at SLAC gives 
the sin Ow with the best precision at low enegries. Fur- 
ther, experimentally M0ller scattering is actively used in 
polarimetry to measure the polarization of the electron 
beams [3| as well as monitoring of luminosity (for exam- 
ple, at DESY Yet another reason stimulating an 
interest in M0ller scattering consists in possibility to test 
the Standard Model and to reveal traces of new physics. 
In the intensively discussed projects of the ILC, e - e~ 
and fj,~ fi~ colliders [5|, high hopes for the discovery of 
Higgs bosons, manifestations of contact interactions, the 
compositeness of the electron, new gauge bosons, etc., are 
pinned on the scattering of identical polarized fermions 
(e,/i). 

A precise comparison of the experimental results with 
the theoretical predictions requires to take into account 
the radiative effects correctly on both QED and elec- 
troweak levels. 

Generally within the polarimetry measurements by the 
M0ller scattering, the value of the transferred momentum 
is rather low, and therefore electroweak effects can usu- 
ally be neglected. But otherwise in the projects of ILC 
where the energies are characterized by the TeV region, 
the value of the weak and electromagnetic effects will 
have the same order. Similarly, due to specific character 
of observables (singly polarized parity- violating asymme- 
try), such experiments as E-158 is sensitive to both elec- 



tromagnetic and electroweak radiative corrections (and 
new physics phenomena at the TeV scales). 

The exact calculation of the lowest-order electromag- 
netic radiative corrections to polarized M0ller scattering 
was performed by Shumeiko and Suarez @. The elec- 
troweak radiative corrections to polarized M0ller scat- 
tering at high energies were computed in Q (without 
hard bremsstrahlung contribution), and at low energies 
corresponding to conditions of E-158 were computed in 
papers [8j (without hard bremsstrahlung) and j9tll0| (in- 
cluding hard bremsstrahlung). The detailed calculation 
presented there demonstrates the significant value of the 
radiative effects that have to be explicitly included both 
in QED and the Electroweak theory predictions. 

In the given paper we considered the lowest order elec- 
tromagnetic radiative corrections both to the longitudi- 
nally polarized cross sections and the doubly-polarized 
parity conserving asymmetry 
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where the first (second) cross section subscripts L and 
R correspond left and right degree of beam (target) po- 
larization respectivelly. Similarly to |(|, we perform our 
calculations within the covariant Bardin- Shumeiko ap- 
proach 0, that allows to cancel out the infrared 
divergences in such a way that the final result does not 
depend on any unphysical parameters (such as a frame- 
dependent cutoff AE that separates the soft photon con- 
tribution region from the hard one). Using the ultrarel- 
ativistic approximation allows us to obtain the compact 
form for radiative correction expression that is conve- 
nient (and sometimes necessary) for fast and more pre- 
cise computer treatment. Moreover during the numer- 
ical estimations it was found that the numerical result 
strictly depends on missing mass cuts. At the same 
time it should be stress that the first correct applica- 
tion of the some kinematical cuts within the covariant 
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FIG. 1: The lowest order graphs giving contribution to the 
M0ller scattering: a) t-channel; b) w-channel. 



Bardin-Shumeiko approach was presented in flU , Il4| . 
Unfortunately Ref. |y] did not investigate the effects of 
experimentally- motivated kinematical cuts. 

The paper is organized as follows. In the Section II the 
kinematics of M0ller scattering as well as the cross sec- 
tion and asymmetry at the lowest order are introduced. 
In the Section III the structure of the lowest order radia- 
tive corrections (virtual and real photon contributions) 
is explained. The Section IV presents numerical results 
applied to the kinematics of E-158 (SLAC) and JLab 
experiments. The Section V contains some conclusions. 
The explicit expressions for finite part of the real photon 
emission could be find in Appendix. 



II. THE LOWEST ORDER CONTRIBUTION 



written as follows 
2na 2 " 



t 2 



(1 + P)--(1-P)- 

s u 



+ (*<-> u). (6) 



Here and later each a denotes the differential cross sec- 
tion over the kinematic variable y (a = da/dy) that is 
defined as 



y 



(7) 



P = PbPt, where Pb,Pt are the polarizations of the 
beam and target electrons. 

The form of the Born cross section (JSJ) with factorized 
combinations 1 ± PbPt is very convenient for construc- 
tion of the polarization asymmetry Q that does not de- 
pend on any energies: 

o _ y{l -y)(y 2 -y + 2) _ sin 2 0(7 + cos 2 B) 

LR - /, i „ /„, T\Y2 _ to i „„„2^2 ' V 8 J 



(l + y(y-l)) 2 



(3 + cos 2 ey 



where 9 is a scattering angle of the detected electron with 
4-momentum k 2 in the center mass system of the initial 
particles k% + pi = 0. The cosine of this angle can be 
express via invariants in the standard way: 



cos6»° = 1 + 2t/s = l-2y, 



(9) 



while the energy of the scattering lepton in Lab. system 
reads: 



The lowest order Feynman graphs giving the contribu- 
tion to M0ller scattering 



e(fci,£t) + e(pi,r)t) -> e(fc 2 ) + e(p 2 ) 



(2) 



are presented in Fig. ^ Here kx, pi (&2, P2) are the 
4-momenta of the incoming (outgoing) electrons (k 2 = 
k 2 = Pi — p\ = m ) while the beam and target (77^) 
polarization vectors read: 
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y/s(s - 4m 2 ) V m 
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k\ - 2mp 1 , 



y/s(s — Am 2 ) 



„ , s - 2m 2 
2mk\ p\ 



(3) 



Then s and other Mandelstam variables can be introduce 
in the standard way: 

s = (ki +pi) 2 , t=(ki-k 2 ) 2 1 u=(k 2 -pi) 2 , 
s + t + u = 4m 2 . (4) 

Notice that for the Born kinematics (strictly speaking for 
the non radiative process) 



u = uq — 4m 2 — s — t. 



(5) 



Neglecting the electron mass, the Born cross section 
for the M0ller scattering of longitudinal electrons can be 
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s+t- 2m 2 
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(10) 



III. ELECTROMAGNETIC RADIATIVE 
CORRECTIONS 

The lowest order radiative corrections to M0ller scat- 
tering appears from the graphs with the additional 
virtual particle (V-contribution, see Fig. [3 for the t- 
channel) and from the real photon bremsstrahlung (R- 
contribution, see Fig. [21 for the t-channcl). It should be 
noted that both these parts include the infrared diver- 
gency but their sum must be infrared free. In this sec- 
tion the explicit expression for V- and R- contributions 
as well as their infrared free sum are presented. 



A. Virtual contribution 

For the calculation of one-loop electromagnetic ra- 
diative corrections we apply the on-shell renormaliza- 
tion scheme of electroweak standard model. The build- 
ing blocks needed for explicit calculations according this 
scheme have been worked out in paper of Bohm et al. 
[l5| . where we take the results for gauge boson self- 
energies and vertex functions. 
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index / = e, fx, r). Hadronic part of the photonic 
vacuum polarization associated with light quarks 
can be directly obtained from the data on process 
e + e~ — * hadrons via dis per sion relations. Here we 
use parameterization of 

ReIL h (-t) S A + Blog(l + C\t\), (13) 






with updated parameters A,B,C in different energy 
regions. 

2. For the contribution of electron vertices we used 
the results of the paper 0] (see also references 
therein). We can obtain the vertex part as 



t 2 



L s 

-(1 - P)— AUt, m 2 ) + (t <-► u), (14) 
u J 
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+ log J-L + log^ + 4(— - 1). (15) 



FIG. 2: The virtual one-loop graphs giving contribution to 
the corrected M0ller scattering within t-channel. 



The virtual contributions to M0ller scattering can be 
separated into three parts: 
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Box 
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3. Our calculation for the box cross section gives com- 
pact formula: 



Box 
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where 

1. a s is a virtual photon self-energy contribution 
[Fig. El (a)]; 

2. <7 1/er is a vertex function contribution [Fig.[2(b,c)]; 

3. <7 So:E is a box contribution [Fig. [21 (d,e)]. 
Now we consider each of them. 

1. The contribution of the virtual photon self ener- 
gies (including the photon vacuum polarization by 
hadrons) to the cross section looks like 



Ana* 



(1+P)- 



(1-P)— + (t«-»u). (12) 



Here £J(— *) is the renormalized transverse part 
of the 7-self-energy (this part includes vacuum 
polarization by e, and r charged leptons: in corre- 
sponding formula of |15| we should take a summing 



; (77) 
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The expressions < have the form 
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and logarithms look like 



log tt, lx = log ■ 
6 |i| 5 t 



(16) 



(17) 



(18) 



It should be noted that vertex and box parts contain 
the infrared divergence through the appearance of the 
fictitious photon mass A. The infrared part from virtual 
cross section can be extracted in a simple way: 



a v (X 2 ^s) 
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It should be noted that due to kinematical restrictions 
the upper limit of the integration in (|21[) is defined as 



st + y/s(s - 4m 2 )t(t - 4m 2 ) 
2^2 



s + 1. (23) 



On the other hand, the energy of the scattering lepton 
in Lab. system I|1U|I for the radiative process transforms 
into 



Et = 



t-v- 2m 2 
2m 



and reaches its minimum value for v = v max 

S + t- V max - 2m 2 



s 2 + e 



2m 



2st 



(24) 



(25) 



FIG. 3: The real one photon emission graphs giving contri- 
bution to the corrected M0ller scattering within f-channel. 



B. Real bremsstrahlung contribution 

The full set of Feynman graphs contributed to the real 
photon bremsstrahlung are presented in Fig. [3J For ex- 
ivergi 

of Bardin and Shumeiko [llj : 



traction of the infrared di verg ency we use the prescription 

(20) 



a R = a R -a 



IR 
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where the infrared free part can be presented in the fol- 
lowing way 



Obviously the electron with the energy ij^STl cannot be 
detected. Moreover, as it was point out first in p|, the 
variable v can be directly reconstructed from the mea- 
sured momenta. However not all events with non-zero 
v < Vcut can be rejected from the experimental data 
due to finite resolution of the experimental equipment. 
Therefore during the radiative corrections calculation for 
the given experimental setup it is necessary to take into 
account this fact. 

Notice that for the radiative events the cosine of the 
scattering angle in the centre mass system of the initial 
particles also depends on integration variable v: 



cosO 1 



1 + 2t/(s - v). 



(26) 
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The explicit expressions for Si are presented in the Ap- 
pendix. The integration in (|21() is performed over vari- 
able v that is a so-called inelasticity. The reason of this 
term can be explain by the fact that for the radiative 
process the last relation in Q transforms into 



s + t + u = v + 4m . 



(22) 



The infrared-divergent part of bremsstrahlung cross 
section integrated over the real photon phase space is 
given in terms of a finite (and infinitesimal) photon mass 
A in 



r « 
a iR = - 

7T 



41og- 



The explicit expression for v can be defined as v = A — 
m 2 , where A = k\ — &2 +Pi and A 2 is a so-called missing 
mass squared. 
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where (see for details) 
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FIG. 4: The relative corrections to the unpolarized (S u ) and polarized (5 P ) parts of the cross section as a functions of the 
scattering angle (defined according to eq. 1331 ) for JLab (Eb ea m = 1 GeV, solid lines) and SLAC (Ebeam = 45 GeV, dashed 
lines) kinematic conditions with different inelasticity cuts: 1) v cut = 0.5v max ; 2) 0.9v max ; 3) 0.99v max . 
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Here Li 2 (:c) is the Spence function and 



l m = log— 5-, l r = log 

TJl 

Summing up ifTTjl and J23 
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we obtain a cancellation of infrared divergencies from R- 
and V- contribution. 

Finally, the total infrared free radiative corrected cross 
section reads: 



on the scattering angle in the center mass system of the 
initial electrons, we investigate the dependence of radia- 
tive corrections on the value of the missing mass cut. 

As it was mentioned above for the radiative events the 
cosine of the scattering angle has to be expressed not 
only via t and s as for non-radiative events iJjJJ but and 
via inelasticity v too (see (|26fl ). Taking into account that 
we calculate the cross section as a function of y or t = 
—ys (because s is fixed) and inelasticity is the inegration 
variable we has some uncertainties in the definition of 
the scattering angle for the observable cross section 13211 . 
In order to escape it we use the standard non-radiative 
approximation i.e. we assume that 

cos 9 R - cos 6° = 1 + 2t/s. (33) 



obs 



a + a 



RV 



(32) 



IV. NUMERICAL ESTIMATIONS 



Basing on the e quat ions presented above the FOR- 
TRAN code MERA[19| (M0ller scattering: Electromag- 
netic RAdiative corrections) has been developed. In this 
section using MERA the numerical estimation of radia- 
tive effects to the M0ller scattering of longitudinally po- 
larized electrons is presented. 

There are two basic differences between the numerical 
analysis that are performed in this and previous 0, ITfj| 
papers: we show the dependence of radiative corrections 



The cross section for polarized M0ller scattering can be 
presented as a difference of the unpolarized and polarized 
parts 

„0,obs „0.obs D^-0,o6s fQA\ 

a = (T u - Fa p > . (34) 

In the Fig. 21 the ^-dependence of the relative radiative 
correction for the unpolarized and polarized parts of the 
cross section 



(35) 

for three different inelasticity cuts v cu t (and, therefore 
A 2 cuts): v cut = 0-5v max , 0.9v ma x and 0.99v ma x is pre- 
sented. One can see the following features of their behav- 
ior: the presence of maximum values at 9 > 90°; sizable 




FIG. 5: ^-dependence (6 is denned according to eq. jSSt ) of the Born (solid line) and observable (dashed lines) asymmetries for 
JLab (Ebeam = 1 GeV) and SLAC (Eteam = 45 GeV) kinematic conditions with different inelasticity cuts: 1) v cu t = 0.5v m ax', 

2) 0.9^ ma xi 3) Q-99Vmax • 
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FIG. 6: The relative corrections to the asymmetry (I36H 
as a functions of the scattering angle (defined according to 
eq. 133H ) for JLab (Ei, eam — 1 GeV, solid lines) and SLAC 
(Ebeam = 45 GeV, dashed lines) kinematic conditions with 
different inelasticity cuts: 1) v cut = 0.5v ma x', 2) 0.9v ma x', 3) 
0.99?; max ■ 



increasing when v cut tending to its maximum value. For 
the v cut = B-5v max the corrections S u (S p ) for 9 = 90° 
equal to 1.075 (1.064) for SLAC and 1.054 (1.045) for 
JLab. 

The (^-dependence of the Born and observable asymme- 
tries with the same inelasticity cuts is presented in Fig.0 
In this figure it can be seen that in the most region of 



9 the corrected asymmetries are less then the Born ones 
and essentially decrease with the increasing v cu t- The 
Fig. El where ^-dependence of the relative corrections to 
the asymmetries 

Aobs _ AO 

§a = A LR A LR (36) 

is presented, reflects this fact clear. Particularly, it could 
be seen that for the realistic v cut — 0.5v max we have 
8 A = -0.008(-0.01) for JLab (SLAC) at 9 = 90° while 
for v cut = 0.99v max the relative corrections to the asym- 
metries reach Sa = — 0.4(— 0.55) for JLab (SLAC) at the 
same angle. 

Let us consider the situation with more realistic small 
v cu t ■ As it could be seen from 129|) and (|31|) the radiative 
corrected cross section (|5"2l diverge when v cu t tends to 
zero. Such cross section behavior can be explain in a 
simple way. Naturally that there is no any real photon 
emission in the limit v cu t —> 0. Therefore we need to 
say about the infrared divergency that appear from V- 
contribution and can not be canceled due to any real 
photon emission absent. 

The other very interesting feature consists in the de- 
viation of the observable asymmetry from the Born one 
at the small v cu t where asymmetry reach its maximum 
value. Due to rather small effects once again in Fig.0the 
quantity i|3b|) as a function of the ratio v cut /v max is pre- 
sented for JLab and SLAC kinematic conditions. From 
this picture it can be seen that for 0.001w maa; < v cu t < 
0.lv max the relative correction to the asymmetry is flat 
and consists some dozen of percent while starting with 
Vcut > 0-lvmax it rapidly falls. 
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FIG. 7: The relative corrections to the asymmetry (I36H as 
a functions of v cu t/v max at the scattering angle (defined ac- 
cording to eq. in CM system 1) = 90°; 2) 100° for 
JLab (Ebeam = 1 GeV, solid lines) and SLAC (Ei, eam = 45 
GeV, dashed lines) kinematic conditions. 

V. CONCLUSIONS 

The explicit expressions for the lowest order electro- 
magnetic radiative corrections to M0ller scattering of lon- 
gitudinally polarized electrons in ultrarelativistic approx- 
imation have been obtained. Basing on these expressions 
the FORTRAN code MERA has been developed. 

The numerical analysis performed for different values 
of missing mass (inelasticity) cut has shown that the ra- 



diative corrections are strongly depended on this param- 
eter. So when this cut tends to its maximum value two 
tendencies should be observed: the radiatively corrected 
cross sections increase while the radiatively corrected 
asymmetries decrease. At the same time ^-dependence 
has shown some common features: the relative correc- 
tions to the cross sections (asymmetries) has a maximum 
(minimum) at 9 > 90°. For more realistic small cuts the 
relative corrections to the asymmetry are rather flat and 
amounts some dozen percent. 

Taking into consideration a large scale of obtained ra- 
diative effects we proof the necessity of radiative correc- 
tion procedure for JLab and SLAC experiments. Partic- 
ularly to perform data processing correctly it is necessary 
to construct Monte Carlo generator for simulation of ra- 
diative events within M0ller scattering. 
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APPENDIX 

Si is the f-channel contribution of the emission from 
the upper electron leg: 



Si = S? + PSf + S%, 

= -L A -L A + L s (2s 2 t- 2 + 2t~ 1 u + 2r 2 u 2 + I) + L x {-2sH~ 2 - 2st _1 - 2r 2 u 2 - I) - 2L t + AL m 

-2t- 2 (v{s 2 + u 2 )s _1 tt _1 - 2t), 
Sf = -L A + L a + L s (-2t~ 1 u - 1) + L x (-2sr 1 - I) + AL m t^ 1 (s - u) + Ar 2 (s + t-u) 

-2(4s 2 + 2st - t 2 )s- 1 t~ 1 (s + ty 1 - 2(3s + 2t)r 1 u~ 1 + 2us~ 1 t~ 1 , 
ST = A(l + P)t- 2 (u + u Q ){l + tL m ); (A.f) 

S2 is the i-channel contribution of the interference from the upper and lower electron legs: 

s 2 = {s^ + ps p 2 )it + s^ 
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S% = L s (8s 3 t _1 - 8s 2 r x v + 8s 2 + Ast + 4si" V - 6sv - t 2 + 2tv - v 2 )/(t - v) + L^SsH' 1 + Ws 2 ^^ 
-16s 2 - 12st - 12si"V + 22sv - 5t 2 + I2tv + Ar x v 3 - Uv 2 )/(t - v) + 2L u (2s + t-v) 
+L 1 (4s 3 t- 1 - As 2 T 1 v + 2s 2 + st + 2st~ 1 v 2 -sv + tv- v 2 ) + L 2 {-AsH~ 1 + SsV 1 *; - 10.s 2 - 9st 
-6si~V + 13sv - 3t 2 + 7tv + 2t~ 1 v 3 - 6v 2 ) + 2L 3 t(-2s - t + v) + 2t((v - s)" 1 - (s + t)" 1 
+2(2s + t-v)/(t- v) 2 ), 

Sf = L s (l2s 2 t + 8s W - 20s 2 v + Ast 2 - 12stv - 4st~ V + I2sv 2 - t 3 + 3t 2 v - 3tv 2 + v 3 )/(t - vf 

+L x (-I2s 2 t - 8s W + 20s 2 v - 20st 2 + h2stv + 12st~ 1 v 3 - AAsv 2 - 7t 3 + 2U 2 v - 33tv 2 - 4f"V 
+19v 3 )/(t - v) 2 + 2L u t + L 1 (2s 2 t 2 - 8s 2 tv - 4s W + 10s V + st 3 - 3st 2 v + hstv 2 + 2st~ 1 v 4 
-5sv 3 + t 3 v - 3t 2 v 2 + 3to 3 - v 4 )/(t - v) 2 + L 2 (-2s 2 t 2 + 8s 2 tv + 4s W - lOsV - 3sf 3 
+17st 2 w - 31sto 2 - 6si" V + 23sv 3 - t A + 7t 3 v - 17t 2 v 2 + I9tv 3 + 2t~ 1 v 5 - 10u 4 )/(i - v) 2 - 2L 3 tv 
-2t(t + v)(s + t)- 1 (s-v)-\ 

5 2 a = 2t- 2 {(L x - L A )(s 2 (l -P) + (u 2 + uu + u 2 ){\ + P)) + s(L A + L s )(l + P)(u + u )}; (A.2) 
53 is the i-channel contribution of the emission from the lower electron leg: 
S 3 = (S^ + PSl)/t 2 + Sl 

5^ = -L u t(s 2 + u 2 ){v - 2t)(t - v)- 2 + {isH 3 ^ 1 + 12s 2 t 2 - 13s 2 to + 4sV + 3st A r- 1 + 9st 3 
-2bst 2 v + I7stv 2 - Asv 3 + ^t' 1 + 1 4 - 13t 3 t; + 19t 2 v 2 - y tv 3 + 2v 4 )/(t - v) 3 , 

Sf = L u t{s-u)(2t-v)(t-v)- 1 + ^—^(st 3 + 8t 2 v-7tv 2 +Av 3 )(t-v)- 2 , 

2t 

SI = 2{l + P)t- 2 (u + u )(l + tL u ); (A.3) 
•$4.5,6,7 are the contributions of the interference between the t- and u-channel graphs: 
S 4 = ((l-P)S r i + PS p 4 )/u + St, 

SI = -2L lS 3 r 1 + L 2 {2s 3 r 1 + 4s 2 + 3si + t 2 ) + 2L A (-'3s 3 - 7s 2 t + 2s 2 v - 5st 2 + 3sto - t 3 

+t 2 v)(s + ty 2 + 2L m (2s - v) + L s sr x (2s - v) + L x {-As 2 r x + 3st _1 w - 6s - M - r x v 2 + 3v) 
+L t (-2s -3t + 2v) + L 3 t(2s + t) + 2{-2s 3 r 1 v + s 3 - 2s 2 t - 3s< 2 + 3stv - sv 2 
-t 2 v + tv 2 )(s + t)~ 2 (t - 

SJ = AsH^vis + t)- 2 , 

St = 2s 2 (l- P)(2-u L A + 2tL m + sL A )/(tuuo); (A.4) 



S 5 = (1-P)SI + S*, 

S% = L 5 v + 2L 6 (~s 3 - s 2 u + s 2 v - sv 2 - uv 2 + v 3 )(s + - uy 1 ^ 1 - L A + 2L m (-2s 4 u~ 1 

+5s 3 u~ 1 v - 4s 3 - 3s 2 u - 6s 2 u~ 1 v 2 + 10s 2 v - 2su 2 + 9suv + 4su~ 1 v 3 - llsv 2 - u 3 + 4u 2 v - 6uv 2 - u~ 1 v i 
+4v 3 )(s + uy 1 ^ - uy 1 ^ 1 + L r (-2s 3 - 2s 2 u + s 2 v - u 2 v)(s + - u) -1 ^ 1 

+L s (4s 4 u~ 1 - 4s 3 u~ 1 v + 4s 3 + 4s 2 u + 3s 2 u~ 1 v 2 - 8s 2 v - 4suv - su~ 1 v 3 + 6sv 2 + 2uv 2 - 2v 3 )(s + uy 1 
(v - uy 1 ^ 1 - L t + L u - L x {4s 3 u + 2s 3 u~ 1 v 2 - 6s 3 v - s 2 u 2 + 2s 2 uv - s 2 v 2 - 2su 3 + 7su 2 v - 8suv 2 + 3sv 3 
-u 4 + 3u 3 v - 4uV + 3uv 3 - v 4 )(s + uy x (u - u)~ 2 i _1 + 4(2s - ^tT 1 * -1 , 
SI = 2s 2 (l- P)(2 + sL s + 2tL m + u L x )t- 1 u- 1 UQ 1 - 1 (A.5) 

S 6 = {l-P)S c J{tu) + S a & , 

S c 6 = 2L u s 2 - 3r _1 s 2 t 3 (t - v)~ 3 + t(5sH - 2s 4 v + 10s 3 t 2 - 4s 3 tv + s 2 t 3 + I0s 2 t 2 v - \2sHv 2 + 4s 2 v 3 - 4st 4 

+Ust 3 v - 18siV + lOstv 3 - 2sv 4 + t 4 v - 3t 3 v 2 + 3t 2 v 3 - tv 4 )(s + ty 2 (t - V y 3 
Sg = 2s 2 (l - P)(sL s - u a L A + tL u )t- 1 u- 1 U Q 1 ] (A.6) 

S 7 = ((l-P)S c 7 + PS^)/(tu) + S^ 

Sy = -2L x s 3 + L^s 2 + t 2 ) + L 5 (s + t)(s -t-v)u + L s s(2s - v) + L x (-4s 2 - 2st + 3sv + tv- v 2 ) 

+L u (-3s 2 t 2 + 4s 2 tv - s 2 v 2 - 2st 3 + 6st 2 v - 6stv 2 + 2sv 3 - t 4 + 3t 3 v - 4t 2 v 2 + 3tv 3 - v 4 )(t - v)~ 2 
+2(-4s 3 t + 2s 3 v - 7s 2 t 2 + 5s 2 tv - 3st 3 + 3st 2 v - t 3 v + 2t 2 v 2 - tv 3 ){s + ty\t - v)~ 2 , 

S'P = 4s- 1 v{s 3 + 2t 2 v-4tv 2 + 2v 3 )(t-vy 2 , 

Sj 1 = 2s 2 (l- P){2 + u L x +tL u + sL A )t- 1 u- 1 u^ 1 ; (A.7) 
•58,9,10 are the pure u-channel contributions: 

5 8 = (S% + PSl)\ t „ u /u 2 + 5 8 Q , 5 8 Q = ~2{s 2 +t 2 - (s 2 - t 2 )P)((u + u )u~ 2 u^ 2 - Lau^Uq 1 ); (A.8) 

S 9 = (S% + PS p 2 )\ t ^ u /u + SI, 5 9 Q - -2(.s 2 + t 2 - (s 2 - t 2 )P)(2tL m + sL A + sL s + tL u )(u + u ) U - 2 u^ 2 - (A.9) 

S 10 = (S? + PSf)\ t ^u + S? Q , S? Q = -2(s 2 + t 2 - (s 2 - t 2 )P)(2(u + u )u- 2 u„ 2 + L^ 1 ^ 1 )- (A.10) 

Evidently, the contributions S\, S2, S3 are in agreement with the corresponding terms of calculations [T^ (un- 
polarized fermion scattering) and 01 (longitudinally polarized fermion scattering), if we suppose that masses of the 
initial fermions being equal. 

Here we present the logarithms (and their combinations) which were used in hard bremsstrahlung calculation (notice 
that all of them do not lead to infrared singularity): 




v — s 



1 



log 




v — u 



1 



log 




v - t 



1 



log 



T{v~t) 2 



m 2 t 2 



10 



1 (v-t) 2 1 

L u = 7 log 1 — , Li = -(L s - L A ), 

v — t m z T v 

L 2 = -{L x + La), L3 = —L t , 
v v 

1 1 u 2 m 2 

L4 = -(L u - 2L m ), L 5 = — — — — log r^—, 

v v(s + t) [v — u) z t 

L 6 = - log L r = - log -3 , t = v + m . 

s m z ut u u z 
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